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Dh; Abstract: We study the Higgs sector of the secluded [/(l)'-extended MSSM (sMSSM) 

r~| \ focusing on CP violation. Using the one-loop effective potential that includes contributions 

from quarks and squarks in the third generation, we search for the allowed region under 
theoretical and experimental constraints. It is found that the possible region for the elec- 
troweak vacuum to exist is quite limited, depending on the parameters in the model. The 
masses and couplings of the Higgs bosons are calculated with/without CP violation. Even 
at the tree level, CP violation is possible by complex soft SUSY breaking masses. Similar 
to the CPX scenario in the MSSM, the scalar-pseudoscalar mixing enables the lightest 



<^ ' Higgs boson mass to become smaller than the Z boson mass while the coupling with the 

00 ' 

o 



Z boson is sufficiently suppressed to avoid the LEP experimental constraints. However, 
unlike the CPX scenario, large fi and A are not required for the realization of large CP 
violation. The typical spectrum of the SUSY particles is thus different. We also investi- 
gate the possible upper bound of the lightest Higgs boson in the case of spontaneous CP 
violation. The maximal value of it can reach above 100 GeV with maximal CP-violating 
phases. 
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1. Introduction 

Many new physics models have been proposed to address the issue of the so-called gauge 
hierarchy problem that cannot be resolved within the framework of the standard model 
(SM). Supersymmetric extensions of the SM have been paid much attention as possible 
solutions to this problem. In particular, the minimal supersymmetric standard model 
(MSSM) can solve not only this problem but also cosmological problems such as dark 
matter and baryon asymmetry of the Universe and so on. Nevertheless, the model still has 
an unattractive feature: the /x problem, where /i appears in the mass term of the higgsinos. 
As long as no special symmetry exist in the theory, the scale of fi is supposed to be the 
grand unified theory (GUT)/Planck scale from the naturalness point of view. However, 
once the electroweak symmetry is broken, the scale of fi should be at about the W boson 
mass. One direction to provide a natural scale for /i is to introduce a gauge singlet field 
{S) into the MSSM. Several variations of this extension have been proposed: the next-to- 
MSSM (NMSSM) |, |, §, the nearly MSSM (nMSSM) @, |, the [7(l)'-extended MSSM 
(UMSSM) I, 0, 1, and the secluded C/(l)'-extended MSSM (sMSSM) [|, 0. Comparisons 
among these singlet-extended MSSM models can be found in Refs. [|Tl|. A common feature 



in these models is that there is no fundamental /i term in the superpotential. After the 
symmetry breaking associated with the singlet field S, the // term is effectively generated 
by the product of the dimensionless coupling and the vacuum expectation value (VEV) of 
S, and thus no fine tuning is required. Because of the introduction of singlet field(s), such 
models have richer physics than the MSSM. 

In this paper, we focus on the Higgs sector of the sMSSM with particular emphasis on 
CP violation. The sMSSM is a string-inspired model whose particle content of the Higgs 
sector comprises two Higgs doublets and four Higgs singlets. They are charged under the 
SU{3)c X SU{2)l X U{1)y x {7(1)q/ gauge symmetry. Once the additional C/(l) symmetry is 
introduced, a new gauge boson Z' must exist in the model and can mix with the ordinary 



Z boson |12, 13|. From the negative results of Z' search at LEP, the magnitude of the 
mixing angle between them (denoted by azz') must be suppressed at 0{W~^) level |14|. 



The sMSSM provides an explanation for such a Z-Z' hierarchy in a natural way. If the 
U{iy symmetry is broken around the TeV scale, the VEVs of the additional three Higgs 
singlets (5i, 52, S3) are expected to be of O(TeV). This makes azz' small enough to escape 
from the current experimental bounds on the Z' boson. 

Due to the extension in the Higgs sector, it is possible to break the CP symmetry 
explicitly and spontaneously at the tree level, which is forbidden in the MSSM. It is well 



known that the Kobayashi-Maskawa CP- violating phase |15| in the SM is too small to gen- 
erate sufficiently large baryon asymmetry of the Universe as observed today [^]. Therefore, 
additional CP-violating phases are required for successful baryogenesis. So far, electroweak 
baryogenesis have been studied in the singlet extended MSSM models: the NMSSM ||17|| , 
the nMSSM [|, |l|], the UMSSM ||] and the sMSSM Q. A detailed analysis of the 



connection between CP violation and baryogenesis, however, is beyond the scope of this 
paper. 

In our analysis, we use the one- loop effective potential that includes contributions 
from the third-generation quarks and squarks. We search for the parameter space allowed 
by imposing both theoretical and experimental constraints on the model. Owing to the 
presence of extra Higgs singlet fields, the tadpole conditions defined by the first derivatives 
of the Higgs potential do not always give the desired vacuum, v = 246 GeV. Therefore, we 
also numerically check whether or not the minimum is located at 246 GeV. We find that 
the possible region for the electroweak vacuum is quite limited, depending on the model 
parameters. 

In the sMSSM, the only source of physical CP violation at the tree level comes from the 
relative phase between the soft SUSY breaking masses and the phases of the Higgs fields. 
We calculate the Higgs boson masses and the couplings between the gauge bosons and 
Higgs bosons in the cases of explicit CP violation (ECPV) and spontaneous CP violation 
(SCPV). It is found that due to the new CP- violating phases, the mass of the lightest Higgs 
boson can be smaller than that of the Z boson. On the other hand, the coupling of the 
lightest Higgs boson to the Z boson is sufficiently suppressed, similar to the CPX scenario 



in the MSSM |21, |2j, |2^. Nonetheless, the fi and A parameters are not necessarily large 
in this model, making the spectrum of SUSY particles different from the CPX scenario. 
We also provide a bound on the lightest Higgs boson mass in the case of SCPV. 



Table 1: Particle content in the Higgs sector of sMSSM 



Higgs 


SU{3)c X SU{2)l X C/(l)y X [7(1)'^, 


Hd 


(i,2,-1/2,Qh/\ 
\ J 


Hu 


{i,2M2,Qh}j 


S 


(1, 1, 0, Qs) 


Si 


(1, 1, 0, Qs,) 


S2 


(1, 1, 0, Qs,) 


S3 


(1, 1, 0, Q53) 



Depending on the mass of charged Higgs bosons, the upper bound can reach above 100 
GeV with maximal CP violation. 

The paper is organized as follows. In Section ^, we introduce the model and define the 
CP-violating phases in a reparametrization invariant way. Theoretical and experimental 
constraints are studied in Section |3|. We examine the effects of CP violation on the Higgs 
boson masses and couplings in Section |^. In particular, the explicit CP-violating case is 



presented in Subsection 4.1 and the spontaneous CP-violating case in Subsection 4.2. The 



discussion about electric dipole moments (EDMs) is presented in Subsection 4.3. Finally, 
we summarize the work in Section |5[ Formulas of the Higgs boson masses are given in 
Appendix |^. 



2. The model 

The particle content in the Higgs sector of sMSSM comprises two Higgs doublets {Hd-, Hy) 
and four Higgs singlets (5, 5i, 52, ^s) Q. As listed in Table [l|, each field is charged under 
the SU{?>)c X SU{2)l x U{1)y x [/(l)g, gauge symmetry. Though it is desirable to have 
U{iy charges (Q's) chosen to make the model anomaly free, a complete analysis of anomaly 
cancellation is beyond the scope of this paper ^ . Neither will we address the gauge coupling 
unification issue here as it requires the knowledge of full particle spectrum in the model. 
Instead, we focus exclusively on the Higgs sector. The model which we are considering is 
extended so that no dimensionful parameter exists in the superpotential W: 



WB 



-dijXSHlHl 



XSS1S2S3 , 



(2.1) 



where A and A5 are the dimensionless couplings. Unlike the NMSSM, the 17(1)' symmetry 
forbids a cubic term in the superpotential which can cause a domain wall problem if the Z3 
symmetry is broken spontaneously. Once the Higgs singlet S develops a VEV, an effective 
/i term is generated by /ieflF = A(S'). Therefore, the scale of fieS is determined by the soft 



^To be anomaly free, exotic chiral supemiultiplets are generally required H, KJ, E3]. For our purpose, 
we assume that they are heavy enough not to affect the phenomenology at the electroweak scale. 



SUSY breaking terms. In Eq. ( |2.l| ) only, there is no interaction between the secluded Higgs 
singlet fields 5*1^2,3 and the two Higgs doublets Hu,d and singlet S. 

The Higgs potential at the tree level is given by the F-, D- and soft SUSY breaking 
terms: 

Vo = Vf + Vd + Fsoft, (2.2) 

where each term reads 

Vf = |Ap{|e,,$j;$i|2 + |5|2($tci>^ + cDt $„)} + \Xs\\\S,S2\^ + \S2Ssf + \SsSi\^), (2.3) 

Vd = ^^(^i<l>rf - $t cl>„)2 + ||cl>tc^„|2 
o Z 

+ i^(^QH^^^^^, + QH^^i^, + Qs\S\^ + Y^Qs^\S,f) , (2.4) 

1=1 
3 

i=l 

-{eijXAxS^'a^i + XsAxsSiS2S3 + m^s^SSi + mls,SS2 + m^^s^S^ + h.c). 

(2.5) 

where g2, gi and (7']^ are the SU{2), U{1) and C/(l)' gauge couplings, respectively. We will 
take g'^ = y/5/3gi as motivated by the gauge unification in the simple GUTs. The soft 
SUSY breaking masses mssi and mss2 are introduced to break the two unwanted global 
C/(l) symmetries. This choice is called Model I, where Qs = —Qsi = —Qs2 = Qs-sf^ and 
Qud + Qhu + Qs = 0- Although the other choice dubbed Model II is also possible, we will 
not pursue it in this paper since there is no room for physical CP-violating phases in the 
tree- level potential |^. The secluded sector {Si, §2,53) can interact with the ordinary ones 
{Hd,Hu,S) through the g[ coupling, mssi and mss2- 

In general, the following five parameters can be complex in the Higgs potential: 

A^A, As^As, mis,, rnls2, ^ks2 ^ ^ (2-6) 

After rephasing the Higgs fields, however, four of them can be made real and only one 
CP-violating phase is physical. In the following, we define the CP-violating phase in a 
reparametrization invariant way. It should be noted that in the UMSSM no physical CP- 
violating phase can survive after rotating the Higgs fields and, therefore, the CP symmetry 
cannot be violated in the tree-level Higgs potential. We parameterize the Higgs fields as 

*, = ,«-. f 7!("^ + '=■< + "•■') ) , ^„ = ,* ( « ) , (2.7) 

\ 4>d J \^ivu + hu+iau) J 

S = —r={vs^hs^ias), Si = ^=-{vs,+hs,+ias,), (i = l-3), (2.8) 



where v = yv"^ -|- t!^ ~ 246 GeV. The nonzero 6's can break the CP symmetry sponta- 
neously. However, the 0's are not independent. Here we define the four gauge invariant 



phases by 

<Pl=^5 + ^5i, ^^2 = ^5 + ^52, '^Z = 0s + 0i + d2, <f4 = Os,+0S2+Os3- (2.9) 

For later convenience, we also define ipi2 = —fi + v'2- The first derivative of the Higgs 
potential with respect to each Higgs field must vanish (tadpole conditions). At the tree 
level, we obtain 

^2 I „2 „, „, I \|2 ^/2 



^/dVo\ 2 , al+gj , 2 2^ n ^uVS ^\M\ 2 ^ 2 ^ ^ s'l n A n 



^ /dVo\ 2 9i+9u 2 2n o^'^^s^lAp 2^ 2-,^9l^ . ^ 

vs\dhs/ ^ ^ ^^' ^ vs ^ ^^^ %5 t's 



— (^^52+^^53) + ^ 

.52\9/.52/" "^ ^^^^ ^52 ^"^"^^^^^^ ^52 '^^^52 



ms, - Reimls^e^ni^ " Re(m|^5,e^^-)^ - i?A,- 

L'52 
+ ^(-|.+-|3) + fQ^2A = 0, 

1 / aVb \ 2 o ^'5il'52 , |As|% 2 , 2 N , 5? ^ . r. 

^\M^/="^^3-^As^;^ + ^(%,+%2) + tQ53A = 0, 
Vu \ Odd I Vd\ aau I 

^^^ = lm{m%s,e'nvs + Im(m|^5^e*'^-)^s, + /as^'5i^53 = 0, 



32 

\dasz 



with 



I\sVs^vS2 = 0> 



A = QhA + QhA + Qsvl + Y, Qs.vh 

i=l 
^X — 7^ ) -'A 



V2 ' ^ V2 



2.10) 
2.11) 

2.12) 

2.13) 

2.14) 
2.15) 
2.16) 
2.17) 
2.18) 
2.19) 
2.20) 

2.21) 
2.22) 



Table 2: Physical Higgs bosons in the sMSSM 





CP-even Higgs bosons 


CP-odd Higgs bosons 


charged Higgs bosons 


CPC 


Hi,H2, H3, H4,H^, Hq 


Ai,A2,A3,A4 


H+,H- 


CPV 


Hi,H2, H3, i/4, Hr,, Hq, Hj, Hs, Hq, Hio 


H+,H- 



Rx, 



V2 



h. 



ImiXsAx^e'^P^) 
V2 



(2.23) 



where (• • •) is defined such that all Higgs fluctuating fields are taken to be zero. Here 
all the Higgs VEVs are assumed to be nonzero. For some parameter sets, however, a 
global minimum can be located at the place where some of the Higgs VEVs are zero. 
Of course, such a minimum cannot be found from Eqs. ( |2.10| )-( p.20| ). We will discuss 
the method of minimum search in Section ^. In the current investigation, we do not 
specify any SUSY breaking scenario. Hence the soft SUSY breaking masses are given by 
the tadpole conditions for the CP-even Higgs fields Eqs. ( p.lO| )-( p.l5| ). After solving the 
tadpole conditions for the CP-odd Higgs fields from Eqs. (|2.16[) -( P^20|) , we find 



h = hs = 0, 



VS2 
vs 



lm{mls^e'''f^] 



-Im(m|^5^e^^-)^. 



(2.24) 
(2.25) 

(2.26) 



The CP-violating phases must satisfy Eqs. ( |2.24 )-( 2.26| ) for the vacuum. As a convention, 
we choose the independent physical CP-violating phase to be ^phys = ^^~s{^SiS2^ ~^ '^12- 

2.1 The mass matrix of the neutral Higgs bosons 

The squared mass matrix of the neutral Higgs bosons is a 12 x 12 symmetric matrix taking 
the form 



li 



H^ A^ 



M^ 



H 



N 



M^ 



N 



Ml 
(Mlp) 



Mlp 
Ml 



(2.27) 



where H^ = {Hq = {hd K hs) hg = (/ig, hs^ hs^)), A^ = (ag = (a^ a„ as) a^ = 
{o-Si 0S2 ^Ss))- The subscripts O and S on h/a denote 'ordinary' and 'secluded', re- 
spectively. In Table g, the physical Higgs bosons in this model are listed for both the 
CP-conserving (CPC) and the CP-violating (CPV) cases. After the symmetry breaking, 
two neutral Nambu-Goldstone bosons G^ and G'^ appear and are absorbed by the Z and 
Z' bosons, respectively. It is straightforward to decouple C^ from the squared mass matrix 
(2.27) analytically by performing the rotation 



cos/? sin/3\ /CO 
— sin P cos P ) \ a 



(2.28) 



1 



We diagonalize the reduced 11 x 11 matrix M% numerically: 



where tan/3 = Vu/vd 

i[ag{rn'^,o,rnl,rn2,rn'l,rn'l,rnl,rnQ,rn'^,rng,rng,mlQ), where rrii < rrii^i [i = 
9) and O is an orthogonal matrix. The explicit expressions for the matrix elements in 



Eq. ( |2.27 ) at the tree level are presented in Appendix]^. 

A complex rrig g and/or a nor 
CP-even and CP-odd Higgs bosons: 



A complex rrig g and/or a nontrivial (pi2 can yield nonzero mixing terms between 



MIp cx Im(m|^5^e 



«'/'12^ 



(2.29) 



This gives rise to broken CP symmetry. A detailed discussion about the CP-violating 



effects on the Higgs masses and couplings will be presented in Subsections ^Jj and ^ 
In the CP-conserving case, 7W|p = and Eq. ( 2.27 ) can be decomposed into two 6x6 
sub-matrices. 

Now we consider the one-loop corrections to the Higgs boson masses. It suffices for the 
current investigation to take into account the contributions of the third-generation quarks 
{t,b) and squarks (^1,2,^1,2)- The one-loop effective potential is given by |26] 



Vi 



Nc_ 
327r2 



E 

q=t,b 




2fhl 



a=l,2 




(2.30) 



which is regularized using the DR-scheme. Here Nc denotes the number of colors, m's are 
the background-field-dependent masses, and M is the renormalization scale. We determine 
M by the condition {Vi) = 0, which implies 



InM^ 



Eo [Ea "ii In ml - 2m^ In ml] 



With the one-loop corrections, the tadpole conditions become 



(2.31) 



dVo 



+ 



N, 



c 



167r2 



E 

q=t,b 




a=l,2 



(2.32) 



where m,"^ = ifh?) and cj) denotes all species of the Higgs fields. The one-loop corrections of 
the third-generation quarks and squarks to the Higgs boson masses have exactly the same 
form as in the NMSSM. The explicit formulas can be found in Ref. Q, 

2.2 The mass matrix of the charged Higgs bosons 

The charged Higgs sector is the same as in the MSSM. Once the /i term in the mass formula 
of the MSSM charged Higgs boson is replaced by the effective // term, //eg = Xvse^^^ /V^, 
we can readily obtain the mass of the charged Higgs bosons in the sMSSM. Its squared 
mass matrix is given by 



{<t>Ut 



Ml 




(2.33) 



At the tree level, it follows from Eq. ( 2.33| ) that 

1 / d-'Vo 



m]j± 



sin j3 cos (3 



m 



w 



+ 



2Rx 

sin 2/3 



vs 



\^\\.2 



(2.34) 



Due to the mixing terms between the Higgs doublets and singlets, the relation between the 



charged Higgs boson mass and the CP-odd Higgs boson mass, 'rnr^± 



m 



w 



+ rriA valid in 



the MSSM, breaks down in general. In the limit of A — > and t;5 — > oo with \vs being 
fixed, frigg^ = m^g^ = and without CP violation, one of the CP-odd Higgs boson masses 
is exactly given by 2R\vs/ sm2f3. The mass relation in the MSSM is recovered in this 
particular case. 

At the one-loop level, the mass formula of the charged Higgs bosons takes the form |22 



rr?H± 



2 , 2i?AW5 



+ 



sin 2/3 

Nc 



2 



167r2 sin /3 cos /3 



/i(m|; 






+ 



mi )(mj 
h(m? ) 

bi 



2 

mf 



b2' 



2mlRtVs 
v^ sin^ /3 



+ 



2mlRbVs 



(m? —mi)(nii — m^ ) v^ cos^ /3 

^ 6i ii 01 12 



/Kl'"^t2 

f{ml m 



Am^ml 2 2^ 



f 2 sin /3 cos /3 " 



(2.35) 



where Rt^b = Re(A^j_;,e*'^^)/\/2, A^b are defined as the trilinear couplings in the soft SUSY 
breaking sector, and f {ml, 1712) is defined by 



f{ml,ml 



mt 



m,o 



rrii 



rriT 



In^ 



?TT,n 



TTln 



In^ 



1 



The explicit form of h{m?) is given in Ref. |27]. As is done in Ref. 
by Eq. (p. 35). Therefore, we take m}j± as an input in our analysis. 



(2.36) 
A\\ is determined 



3. Allowed region 

Finding an acceptable minimum of the Higgs potential is a nontrivial task even at the tree 
level. Even if we require the tadpole conditions and positive-definiteness of the squared 
masses of the Higgs bosons, the global minimum can be found at f 7^ 246 GeV. This 
is because of the presence of the Higgs singlets in the Higgs potential. In Ref. Q, the 
following method is adopted to search for the electroweak vacuum. First, the soft SUSY 
breaking masses and the two trilinear A terms {Ax and Ax^) are taken at arbitrary values. 
After finding a viable minimum, all the given dimensionful parameters are rescaled so 
that V = 246 GeV. In this method, all the Higgs VEVs are determined through the six 
tadpole conditions ( 2.1C| )-( |2.15 ). Therefore unlike the MSSM, tan/3 is an output. Our 
method is equivalent to that, but the other way around. Explicitly, we take the Higgs 
VEVs as the inputs, and then perform the minimum search. That is, v = 246 GeV is 



given in advance. However, as we will see in what follows, the desired electroweak vacuum 
does not always exist. For some input parameters, the location of v = 246 GeV can 
be unstable and the true minimum would roll down to another point that does not give 
V = 246 GeV. Redefining such a minimum as f = 246 GeV by rescaling the Higgs VEVs is 
then inconsistent with the original value of tan/3 that is scale independent. Therefore, we 
discard such cases and keep tan/3 as a fixed input. Before showing the numerical results of 
the minimum search, we consider theoretical and experimental constraints in the following 
two subsections, respectively. 

3.1 Theoretical constraints 

The effective potential at the tree level is 

/T/\ ^ 22, 22, 22, Y^ ^22 

i 

o 9 

-RxVdVuVs - RxsVSiVs^vss + ^3^ ^ {vj - vlf 



1^ /'„.2„,2 , „,2„,2 I „,2„,2^ , \^S\ (^^2 „,2 , „,2 „,2 , „,2 „,2 ^ , 9l 
4 



+^(^d< + v^s + viv's) + ^K^^i + v'sA, + ^sAr) + ^^ • (3-1 



In each direction of vs = vs^ and vs = vs2 with other VEVs being zero, we demand 
the coefficients of the quadratic terms be positive so that the effective potential is not 
unbounded from below: 

m| + m|-2Re(m|s,e''^")>0, i = l,2. (3.2) 

Next we consider the vacuum of the Higgs potential. From the tadpole conditions Eqs. (|2.10| )- 



(2.20), the vacuum of the tree-level potential takes the form 



(Vb)vac = -R\VdVuVS + -RxsVs.VS^VSs - ^ ^^ \ vj - vlf 

-^ / 2 2 , 2 2, 2 2\ \^S\ / 2 2 , 2 2 , 2 2 \ 5l a 2 /o o^ 

After eliminating R\ with Eq. ( p.34|) and imposing (Vo)vac < 0, the upper bound on the 
charged Higgs boson mass is obtained: 

2 2 2\\\'^vl 2 .2^0 4i?As 

^H± <mw+ . 2na + ^Z cot 2/3 - . 2 VS^Vs^VSa 

sm 2p v^ sm Ip 

As an example, we plot the maximal value of the charged Higgs boson mass as a 
function of R\g in Fig. 0. We take A = —0.8, \s = 0.1, ^5 = 300 GeV, vs^ = vs2 = ^Sa = 
3000 GeV, and tan/3 = 1 (red solid line), 5 (green dotted line), 10 (blue dashed line). The 
CP-violating phases are assumed to be zero. Since the dominant terms are proportional 
to 1/sin^ 2/3 in m^±^, tan/3 = 1 gives the smallest m^±^ for a fixed Rx^. For Rxg > 0, the 
value of m^'± decreases as Rxg increases. We find a maximum of Rx^ — 640 GeV. 



50 
45 
40 
35 



"^ tan/5 = 10 




100 200 300 400 500 600 700 

Rxs (GcV) 

Figure 1: The maximum of cliarged Higgs boson mass as a function of R\g. We take vs = 300 
GeV, vsi — vs2 = vs3 ~ 3000 GeV, and tan/3 = 1 (red solid line), 5 (green dotted line), 10 (blue 
dashed line). 



3.2 Experimental constraints 

The U{iy charges of the Higgs fields can be constrained by the experimental results of the 
Z' boson search, namely, the lower bound on the Z' boson mass and the upper bound on 
the mixing angle between the Z and Z' bosons. The squared mass matrix of the Z and Z' 
bosons takes the form 



M 



nir 



ZZ' 



•mzg'i {QHa cos^ P - Qhu sin^ P)v 



mzg'iiQHa cos^ 13 - Qh^ sin^ (3)v 



rrir. 



where 



2 92+91 2 

mz = z V , 



ml = 9? (Ql^d + Q\A + Qlvl + E ^''sA 



(3.5) 

(3.6) 
(3.7) 



The eigenvalues of the squared mass matrix and the mixing angle between the Z and Z' 
bosons are respectively given by 



rrir, 



z + ml + V ("^i - "^i')' + 9xk9i + 9{){QhA - Qhu^D 



.,2\2 



a 



ZZ' 



arctan 



'^mzg'iiQHd cos^ f3 - Qh^ siv? p)v 



rrir 



rrir 



(3.8) 
(3.9) 



The experimental constraints on the Z' boson are rather model-dependent. Here we 
adopt the typical bounds, m^/ > 600 GeV and azz' < 0(10"^) @. In Figs. |, we plot 
the mz' = 600 GeV contour and curves for azz' = i)--, 3, 5) x 10~^ in the Qhu-Qh^ plane. 
The other U{1)' charges are determined by the gauge invariance and the condition for 



10 



breaking the two unwanted global U{1) symmetries as discussed above. Here we show two 
examples: (A) vs = 300 GeV, vsi = vs2 = vss = 3000 GeV with tan/3 = 1 (upper left 
figure) and tan/3 = 50 (upper right figure); (B) vs = 500 GeV, vs-^ = vs^ = 100 GeV, 
vs2 = 3000 GeV with tan/3 = 1 (lower left figure) and tan/3 = 10 (lower right figure). The 
red dotted lines give the mz' = 600 GeV contour, and the region in between represents 
mz' < 600 GeV. The figures also show curves for azz' = 1 x 10^'^ (dashed line in green), 
azz' = 3 X 10"^ (dotted line in blue) and azz' = 5 x 10^'^ (solid line in magenta). In the 
region where Qh^ ^^^ Qhu have the same sign, the two terms in the off-diagonal elements 
of ■M'^z' fsnd to cancel with each other. The upper right figures show that the tan/? 
dependence on Z' search constraints is rather mild since the denominator in Eq. (|3.9| ) is 
relatively large for case (A). In the lower left figure, the covered areas of quadrants II and 
IX have azz' > 1 x 10^^. On the other hand, large portions of quadrants I and III are 
not strongly constrained. If we take tan/3 = 10, the contours of azz' is distorted and the 
region around Qh^ — Qh^/ ^SiV? (3 becomes allowed. In our numerical study, as long as 
one of vs^ (i = 1 — 3) is taken to be at the TeV scale and Qu^i — ~Qh^ does not hold, 
the constraints from the Z' boson search can be easily avoided. This supports the original 
motivation for the sMSSM as mentioned in the Introduction. 

According to the LEP experiments, the mass of the SM Higgs boson should be larger 
than 114.4 GeV at 95 % CL [0]. However, this lower bound cannot be directly applied 
to models beyond the SM due to the modification of the Higgs coupling to the Z boson 
idnzz)- When the Higgs boson masses are smaller than 114.4 GeV, we require instead 

e < HmnJ , (3.10) 

where £, = 9hzz/9hzz ^^'^ ^ ^^ ^^^ ^^ *^ ^^ upper limit on the HZZ coupling and a 



function of the Higgs boson mass |28, ^. In our analysis, we do not consider the processes 
e^e'^ -^ Z* ^ HiHj. They are expected to be less severe in comparison with the processes 
e+e- ^ Z* ^ H,Z. 

We also consider the Z boson decays, Z — > HiHj and Z — > Hil^l~ for the light Higgs 
bosons, and require that: 

^r(Z ^ HiHj) + Y^T{Z^ H,l+l-) < AFz , (3.11) 

where ATz = 2.0 MeV is the 95 % CL upper bound on the possible additional decay width 
of the Z boson Q. 

The other experimental constraints come from the lower bounds of the SUSY particles. 
The mass matrix of the charginos has the same form as in the MSSM if we replace fi with 

_ / Ma -V2mwcosp\ 

where M2 is the SU{2) gaugino mass. The physical CP-violating phase is 9m2 + ^a + ^'i-, 
where 9m2 ^-nd Ox denote the arguments of M2 and A, respectively. For the lower bound 
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Figure 2: Constraints from the lower bound on mz' and the upper bound on azz' in the Qh^- 
Qud plane. We take vs = 300 GeV, vs^ = fSa = "^Sa = 3000 GeV with tan/3 = 1 (upper left) 
and tan/3 = 50 (upper right), and vs = 500 GeV, vs^ = wsa = 100 GeV, vs^ = 3000 GeV with 
tan/3 = 1 (lower left) and tan/3 — 10 (lower right). 



on the lightest chargino mass Xi i we require rn~± > y/s/2 ~ 104 GeV, where y/s is the 



Xi 



center-of-mass energy at LEP2 [31|. On the other hand, the mass bound on the neutralino, 

-6 
In 



niyO > 46 GeV given in Ref. |14| is rather model-dependent. In fact, it is found that m^o 
GeV is allowed in the i?-parity conserving MSSM without gaugino mass unification 



the sMSSM, the lightest neutralino can even be massless, almost a singlino |33]. Therefore 
we will not put an explicit lower bound on the mass of the lightest neutralino, and not 
require that the lightest neutralino be a candidate for the cold dark matter of the Universe 
as well. 

Now we consider extra contributions to the p parameter. It can be easily shown that 
if a model has only Higgs doublets and singlets, /? = 1 at the tree level. As discussed 
before, as long as azz' < 0(10"^), the deviation of the p parameter from unity due to the 
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Z' boson is small enough to evade the current experimental bound Ap < 2.0 x 10^^ |14|. 
Let us consider the one- loop corrections, focusing particularly on the contributions of the 
physical Higgs bosons rather than including all SUSY particles. The correction to the p 
parameter is given by 



Ap: 



nL(o) 

2 

rrir 



2 

mi 



(3.13) 



where 11^^(0) iy = Z, W) are the transverse parts of the weak boson self-energies at the 
zero momentum. The Higgs boson contributions at the one-loop level take the form 



Ap 



Higgs 



Gp 

8^2^ 



^9H,H,zBb{mH,,mH^) - ^ \9H,Hw\'^B5{mH±,mH,] 



i<j 



,(3.14) 



with 



B^{mi,m2) = 


1,2 2\ m^lml , m? 

1 „(?Tli+m2 + 2 2^'^ 2 

= < 2 mj — TTig m,2 




(mi = m2) 


9h^HjZ = 


-- (OiiOrj - OijOn) sin/3 - i02i07i - 



(mi / m2), 



02jOii) cos 13 



QHiHW 



02i COS 13 - OiiSinfi -iOri, 



(3.15) 

(3.16) 
(3.17) 



where Gp = l/iV^v"^) ~ 1.166 x IQ-^ (GeV)"^. Unlike the MSSM, the custodial SU{2) 
symmetry does not guarantee Ap ^^^^ = due to the contributions from the Higgs singlets. 
Finally we comment in passing on the constraints from B physics. The experimental 
results of Eg -^ fi^fi^, b ^ s^ and B^ -^ T~Vr can give a significant restriction on the 
parameter space. However, so long as we limit our interest to the low tan/3 region (< 20), 
constraints from the branching ratios of Bg — > IJ'~^p~ and B^ -^ tVt- are less stringent. The 
b ^ s'y process can be important for the light charged Higgs bosons scenario, m.jj± < 300 
GeV, in which case the contributions from the charged Higgs bosons and those of the 



charginos have to cancel |34] in a way to be consistent with the data |35]. We leave the 
detailed analysis to another paper. 

3.3 Numerical evaluation 

Now we show the numerical results of the allowed regions in both case I and case H. We 
take 



Qh, = Qh^ = 1, A 



As 



mg = 1000 GeV, 



m. 



tR 



A\ {mH± ) , At = Ab = fies/ tan /3, 
Ma = 200 GeV, 



mr = 500 GeV 

Or 



(3.18) 



where mg, m^ and m^ are the soft SUSY breaking masses of squarks. It should be noted 
that Ax is a function of m,jj± , as given by Eq. ( p. 35 ). In Fig. |3|, the allowed region is plotted 
in the As'-A plane (left figure) and tan/3-mj:^± plane (right figure). The input parameters 
in Case I are 



Case I : m,gg_^ 



^ls2 



(500 GeV)^ m\s.^ = -(50 GeV)^ 



vs = 300 GeV, vs-^ = vs2 



vs. 



3000 GeV. 



(3.19) 
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Figure 3: The allowed region in the A5-A plane (left figure) and tan (3-mjj± plane (right figure). 



We take Qh^ = Qh^ — 1, ^ssi — "*SS2 ~ i^^^ GeV) , m;^^ 



S2 



-(50 GeV)^ vs = 300 GeV, 



«Si = "S2 = 1^53 



3000 GeV. 



For the moment, all the CP-violating phases are assumed to be zero. In the left figure, 
we take tan/? = 1 and mj{± = 300 GeV. All the Higgs boson masses are non-negative in 
the region between the two blue curves. For fixed A, the depth of the vacuum decreases 
as Xs decreases and eventually becomes higher than the origin, as can been seen from 
Eq. (|3.3| ). The dotted curve in magenta corresponds to the critical situation, below which 
the vacuum becomes metastable. The region to the right of the dotted-dashed line in green 
has been excluded by the condition ( p.lOl ). Likewise, the region to the right of the dashed 
line in red is excluded by the chargino lower mass bound. In the right figure, we take 
A = —0.8, Xs = 0.1. As in the left figure, m'jj > is fulfilled between the two blue curves, 
within which the vacuum becomes metastable below the dotted curve in magenta. The 
region below the dotted-dashed curve in green is excluded by the condition ( |3.10| ), and 
that below the dashed curve in black by Ap > 2.0 x 10^'^. Since the Higgs singlets can 
affect the lightest Higgs boson mass, the possibility tan/3 = 1 excluded in the MSSM is 
experimentally allowed in our model. On the contrary, the allowed region is much more 
restricted by the conditions for the desired electroweak vacuum. 
In Fig. ^, we consider 



Case II : mlg^ = (306 GeV)^ m^s^ 
Vs = 500 GeV, vsi = vs^ = 



= (56 GeV)2, ml^s2 = (100 GeV)^ 
100 GeV, vs^ = 3000 GeV. 



(3.20) 



In the left figure, we use tan/? = 1 and mjj± = 600 GeV. The region to the left of the 
blue line is excluded by m'jj < 0, and that above the dashed curve in blue results in the 
situation where ^ = Vq + ^1 is unbounded from below. In the region between the two 
lines in magenta, the vacuum is correctly located at t; = 246 GeV. However, the region 



to the left of the dotted-dashed line in green is excluded by Eq. (3.10). The fact that 
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Figure 4: The allowed region in the As- A plane (left figure) and tan /S-tti// ± plane (right figure). 



We take Qh^ = Qh^ = 1, mlg^ = (306 GeV)^, ml^^ 
GeV, vs, = -^52 = 100 GeV and vs., = 3000 GeV. 



(56 GeV)2 



m 



S1S2 



(100 GeV)2, vs = 500 



mff± in this case is larger than Case I implies that Rx is larger. A small A can make the 
vacuum metastable, as can be seen from Eq. (^). In the right figure, we take A = 0.8 and 
A5 = 0.1. The allowed region is inside the two dotted-dashed curves in green and the two 
dashed hues in orange. The dotted-dashed curves in green are obtained from the critical 
value of the LEP bound (|3.10| ) explained above. The dashed lines in orange correspond to 
azz' = 1 X 10~^. The parameter space is highly constrained in Case II. 



4. CP violation 

In this section, we study the effects of CP violation in the Higgs sector. In the MSSM, the 
CP-violating phase in the Higgs potential can be rotated away by a field redefinition. Hence 
there is no explicit CP violation at the tree level. However, once the one-loop corrections 
from the squark sector to the Higgs boson masses are taken into account, mixing terms 
between the CP-even and CP-odd Higgs bosons are generated. In a specific CP-violating 
case called the CPX scenario, the effects of CP violation is extremely enhanced, and the 



Higgs phenomenology is drastically changed |21, 22, 23]. The lightest Higgs boson mass, 
for example, can become much smaller than the current LEP lower bound due to the large 
■M'gp in the squared mass matrix. Its coupling to the Z boson, however, can be sufficiently 
suppressed to escape from the LEP constraints ||2^. Studies of ECPV have been done in 
the NMSSM I, 



1 13], nMSSM [§ and the UMSSM ||, |3|] as weh. Here we discuss both 
ECPV and SCPV in the sMSSM. 



4.1 Explicit CP violation 

As discussed in Section g, there is one CP-violating phase that cannot be removed by 
rephasing the Higgs fields. In fact, the nonzero CP-violating phases are related to each 
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other in the vacuum through the tadpole conditions for the C-P-odd Higgs fields. At the 
one-loop level, we find 



h 

hs 



Nc 






0, 

Im(m|^52e 



Im(m2 e^^^) 



SS2^^"V = -Im("^liS2e 



VS 



(4.1) 

(4.2) 
(4.3) 

(4.4) 



where Itfi = Im^XAtfic'^'^^) / V2. If It or /^ is nonzero, I\ can be nonzero as well at the one- 
loop level. Nevertheless, we will focus exclusively on CP violation peculiar to the sMSSM, 
and take /j = /f, = in what follows. Since we have the relation Eq. ( |2.35| ), the sign of Rx 
is determined through 

I \|2 



sgn{Rx 



sgn mjj± 






Am 



H± ' 



(4.5) 



where Am^± denotes the one-loop correction to the charged Higgs boson mass. On the 
contrary, there is a sign ambiguity in R\g at this stage. The positivity of the squared mass 
of the Higgs bosons gives us R\g > in most of the parameter space. Now let us define 
Ossi = Argimls^), ess2 = Arg(m|5j, 63,82 = Arg(m|^5j. FVom Eqs. (|4j) and (U), it 
follows that 



dssi 



&SS2 



sm 



sm 



m 



S1S2 
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SSi 



VS2 



VS 



sin(6's,52 + ifu) 



V^i, 



m 



S1S2 



m 



SS2 



1)q 

— - sin(6'si52 + fu) 
VS 



^2- 



(4.6) 



(4.7) 



It should be noted that the arguments in the arcsines should be smaller than one, imposing 
additional constraints on our input parameters. 

The CP-violating phases show up in the mixing terms between CP-even and CP-odd 
parts in the squared mass matrix ( |2.27 ). Let us parameterize ■M.'sp in terms of 3 x 3 block 
entries: 
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hi h^s MIp 
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(O) 
SP 
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SP 
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SP 



(4.8) 



After the conditions (|4.3D and (|4.4| ) are applied, the entries are 



M 



(O) 
SP 



03x3, M 



(OS) 
SP 



lm{ml^g^e 
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(4.9) 



(4.10) 
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Figure 5: The effects of the CP-violating phase on tuh and g^yy. We take mjy± 



tan/3 = 1, 



'SSil 



(306 GeV)2 



'■SS2 



(56 GeV)2, 



'SiSal 



= (100 GeV)2, vs 



600 GeV, 
500 GeV, 



100 GeV, and vs^ = 3000 GeV. 



If ■M.'sp has a large portion in M'j^, the CP-violating effects on the Higgs boson masses 
can be enhanced. To achieve this, we assume large values for Im(r7T,| ^ e^'^^'^)vs2/'vs and 
Im(m|^^^e*'^i2)t;5^/ws' under the conditions (4^) and (|4.7|) , rendering 



\m 



SSi 



\m 



SS2\ 



I 2 I ^52 
I 2 I f^5i 



(4.11) 
(4.12) 



for sin(0s'^52 +9^12) — 1- For the moment, we only consider ECPV, and hence V'l = ^^2 = 0. 



We present two examples: one being Case II as given in Eq. ( 3.20 ) and the other being 
Case III specified by 



Case III : mlg, = (72 GeV)^ m^s^ 
Vs = 300 GeV, vsi = vs^ 



- (280 GeV)2 
1500 GeV, VS2 



"^llS2 



100 GeV)2, 
100 GeV. (4.13) 



We take tan/3 = 1 and mjj± = 600 GeV for Case II and tan/3 = 1 and mj{± = 300 
GeV for Case III. In Fig. ^, we plot mn^ and gjj.yy (i = 1 — 3) as functions of 6siS2 ™ 
Case II. In the CP-conserving case, 0siS2 — 0) the second lightest Higgs boson is CP-odd 
because gn^vv ^^ zero. Around Qs\S2 — ^0°, ii\ and Hi switch with each other and their 
CP characters are exchanged, as can be seen from the right figure in Fig. |5|. As in the 
CP-violating MSSM, due to the large off-diagonal terms A^^p, Hi can become lighter than 
114.4 GeV for ^5^52 ^ 60° with g\ yy being highly suppressed. This possibility cannot 
be excluded by the LEP experimental results. This does not seem to be typical in the 
CP-violating NMSSM [^. Although all the Higgs boson masses are positive in the range 
93° < ^5i52 % 102°, the vacuum is metastable and is thus excluded. In Fig. |6|, we plot 



niH, and g^^ 



vv 



1 — 3) as functions of 6*5152 foi' Case III. When ^5^52 = 0, Hi is the 
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Figure 6: 

tan/3 = 1, 
VSi = Vs., -- 



The effects of the CP- violating phase on rriH and g^vv- We take mH± — 300 GeV, 
\mls,\ = (72 GeV)^ |m|sj = (280 GeV)^, Iml^^J = (100 GeV)^, vg = 300 GeV, 
-- 1500 GeV, and vs^ = 100 GeV. 



CP-odd Higgs boson since (7^ yy = 0. In this parameter set, H3 is the SM-hke Higgs 
boson, corresponding to the decouphng Hmit in the MSSM. Both Hi and H2 are composed 
of almost singlet components. The mass rriHi is always smaller than the LEP bound when 
we vary Os^Si^ ^^^ can become as low as 20 GeV around 03^82 = 102°. Since gjj yy is 
less than 10~^, the associated production cross section of Hi with gauge bosons is highly 
suppressed. The masses and couplings of the other Higgs bosons are not much affected by 
CP violation. 

4.2 Spontaneous CP violation 

In this subsection, we discuss the SCPV scenario. If the model contains two Higgs doublets, 
one of the Higgs VEVs can be complex in principle. In the MSSM, there is no room for 
the relative phase between the two Higgs doublets in the potential in the SUSY limit due 
to C/(l)pQ. The only place where the relative phase can show up is the quadratic mixing 
term between the two Higgs doublets to break the SUSY softly. After imposing the tadpole 
conditions, such a phase disappears. It is found that the one-loop corrections to the Higgs 
potential can induce radiative SCPV [^. However, it leads to the appearance of a light 
pseudoscalar [niA < 6 GeV), which is already excluded by the LEP experiments. Many 



studies have already been done for SCPV in the NMSSM with a Z3 symmetry |41, 42 
43, 44]. According to Romao's No-Go theorem [^, with certain radiative corrections in 
the Higgs sector the condition for SCPV leads to a negative squared-mass mode in the 



Higgs spectrum. However, it is pointed out by Babu and Barr |43] that the large radiative 
corrections from the top/stop loops have not been taken into account in the proof of the 
No-Go theorem. The original saddle point in the Higgs potential can become a minimum 
in this case and, therefore, the tachyonic mode no longer appears. In Ref. [^], the upper 
bound on the lightest Higgs boson mass is found to be about 140 GeV in the case of SCPV 
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Figure 7: The representative solution for non-zero ipi and (p2- 

where the full one-loop corrections of top/stop have been included in their calculations. 
In the NMSSM without a Z3 symmetry, the No-Go theorem cannot be applied any more. 
Hence, the SCPV scenario is viable even at the tree level 



In the sMSSM, SCPV is induced by the nonzero 0's that appear in the quadratic terms 
of the Higgs potential. This is also free from the No-Go theorem. To simplify our study, we 
assume that fn^g^, rn'^g^, "^5152' ^^^ ^^'^ ^S^Xs ^^^ ^11 ^^al. Prom the tadpole conditions 
Q-Q, we find 

a sin (/9i + 6 sin c/?2 = 0, (4-14) 

a cos 991 + & cos C/92 = , (4.15) 

c 

^3 = ^i = 0, (4.16) 



where a = iri'gg^vsvsi, b = ''tT'Isj ^•S'^'S'a ' ^^*^ '^ ~ "^li 52^-^1 '^•^2- When Eqs. ( 4.14| ) and ( 4.15| ) 



have solutions, they form a triangle as depicted in Pig. |7[ The analytic solutions can be 
easily obtained: 

1 f be C b\ /, -,r,x 

cos^, = -^^----j, (4.17) 

1 /ac a c\ , . . 

cos^2=2(^----)> (4.18) 

, , 1 [ ab b a\ ,^ ^^, 

cos(c/5i -^^) = -\ — -, (4.19) 



which give the CP-violating extremum. The Higgs potential has the CP-violating mini- 
mum when acjb < 0. 

We can set 61 = 63^ = without loss of generality in Eq. (|2.9| ). Since (^3 = 994 = 0, it 
follows that 

02 = --{ipi+ip2), es = -iipi+ip2), (4.20) 

^Si = ^('/'i - V2), Os, = --{^1 - ^2)- (4.21) 

We examine the possible maximal value of niH in the case of SCPV. Since the numerical 
minimum search is rather time-consuming, we do not conduct a complete parameter scan. 
Instead, we restrict ourselves to scan only the three soft SUSY breaking masses in the 
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Figure 8: The left plot shows the upper bounds on the four light neutral Higgs boson masses, 
m^^^ (cross in red), rrf^^ (triangle in green), mf^^^ (circle in blue) and m^^'^ (square in yellow), 
as functions of mff± . The right plot shows | sin ipi \ and | sin ip2 \ ■ The crosses in red are for | sin ipi\, 
and the triangles in green for | sin (p2 \ ■ 



following ranges: 






"T-ls: 



(10 GeV)2 - (1000 GeV)^ 
1000 GeV)2 - (10 GeV)^ 



(4.22) 



for fixed values of mfj±. The remaining parameters are chosen as A = —0.8, \s = 0.1, 
tan/9 = 1, Vg = 300 GeV, and Vg^ = Vg^ = Vg = 3000 GeV. In Fig. ^, the maximal 
values of m/^. (i = 1 — 4) (left figure) and | sin 931 1 and | sin 932 1 (right figure) are plotted 
as functions of mjj±. For each fixed mjj± 
(m^ 



SSi ' '"^552 ' 



,m 



all m^^^ are obtained for different sets of 
g s)- One can see that the upper bounds on m//. strongly depend on mfj± 
except for rriH^ . It is found that the upper bound on the lightest neutral Higgs boson mass 
j^max ^g below 125 GeV and can reach up to around 123 GeV for mH± = 334 GeV. Since 
the lightest state Hi is mainly composed of the singlet states, mH^ do not increase even 
if we change the values of {rriq^mi^^mi )=(1000, 500, 500) GeV into, say (3000, 1500, 
1500) GeV. In this case, the second lightest Higgs boson H2 receives corrections from the 
top/stop loops. From the right figure of Fig. ^ one can see that the CP symmetry is 
maximally violated when m^^^ > 100 GeV. 

It is noticed that the CP-violating solutions ipi and (p2 are obtained by solving the 
necessary conditions for SCPV, Eqs. (4.14) and ( 4.15| ). In order to check whether they 
give CP violation at the vacuum, we perform the minimization in the ten-dimensional 
parameter space {vd, Vu, vs, vsi, vs2, vs^, 6*2, 6*5, 6si, 032) y and find that the solutions 
obtained above indeed give the CP-violating vacuum. 

4.3 EDM constraints 

The CP-violating phases can also be constrained by the upper bounds on electric dipole 
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moments (EDMs) of electron, neutron and mercury atom |46, ^. Similar to the MSSM, the 
SUSY particles- mediated one-loop diagrams contribute to the EDMs. However, we assume 
that the only sources of CP violation come from 9siS2 ^^ ECPV and ipi [i = 1,2) for 
SCPV in the sMSSM. Therefore, their contributions to the EDMs generally vanish. At the 
two-loop level, however, the Higgs bosons with indefinite CP properties can contribute to 
the so-called Barr-Zee type diagrams [^] and become sizable when tan/? is large. Since we 
take tan /3 = 1 in the CP-violating cases, we expect that they do not put severe constraints 
on 6*5^ 5-2 or Lpi {i = 1,2). 



5. Conclusions 

We have studied the Higgs sector of the sMSSM with particular focus on CP violation. 
The masses and couplings of the Higgs bosons are calculated using the one-loop effective 
potential, including corrections due to the third-generation quarks and squarks. Imposing 
both the theoretical and experimental constraints, the allowed region is obtained for Case 
I and Case H defined in the text. In short, all Higgs VEVs of the secluded Higgs singlets in 
Case I are taken to be of ©(TeV), and in Case II two of them are of 0(100 GeV) and the 
other of C(TeV). Due to the corrections from the Higgs singlets, the tan/3 = 1 case cannot 
be ruled out by the LEP experimental results. However, the conditions for the desired 
electroweak vacuum generally render a very restrictive parameter space. 

In this model, ECPV can be induced by the nonzero phase of m-|. ^ at the tree level. 
It is found that a large value of Os^S2 can make the lightest Higgs boson lighter than the 
LEP bound of 114.4 GeV, provided that the Higgs coupling to the Z boson is sufficiently 
suppressed, similar to the CPX scenario in the MSSM. Nevertheless, large /i and A terms 
are not required in the sMSSM for the realization of large CP violation. Therefore, the 
spectrum of SUSY particles is generally different from the MSSM CPX scenario. 

We have also investigated the SCPV scenario. Unlike the MSSM, SCPV can occur at 
the tree level in the presence of the nonzero 0's residing in the quadratic terms of the Higgs 
potential. Our analysis shows that in this case the lightest Higgs boson mass has a certain 
upper bound, depending on the charged Higgs boson mass. In a specific case, the maximal 
value of niH-i ^^ around 125 GeV for mij± = 334 GeV with the CP-violating phases being 
nearly maximal. 

In this paper, it is assumed that the only sources of CP violation come from the 
Higgs sector. Such CP-violating phases show up in the Higgs boson-mediated two-loop 
diagrams that contribute to the EDMs of electron, neutron and mercury atom. However, 
these diagrams are not important as long as tan/3 = 1. 



As pointed out in Ref. [20|, a strong first order electroweak phase transition is possible 
in the sMSSM due to the presence of the trilinear term \A\S^d^u- In this case, the light 
stop is not necessarily lighter than the top quark as required in the MSSM. A devoted 
study of the electroweak phase transition with/without CP violation will be presented 
elsewhere 
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A. The mass matrix of the neutral Higgs bosons at the tree level 

Here we present explicitly the tree-level squared mass matrix elements for the neutral Higgs 
bosons. The CP-even part is given by 
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The CP-odd part is given by 
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The mixing between CP-even and CP-odd parts is already given in the main text. 
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